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\ Abstract. We study extractions of |V„ S | based on finite energy sum rule (FESR) analyses of 

hadronic T decay data. We show (i) that the "(0,0) spectral weight" implementation (proposed pre- 
viously in the literature as a favorable version of this analysis) suffers from significant convergence 
K-j , problems, but (ii) that alternate implementations exist which bring these problems under control. 

■ Results based on present spectral data are shown to be in agreement with those of other approaches, 

though with somewhat larger experimental errors. Sub- 1% determinations of |V„ S | are also shown to 
be possible once T data from the B factories becomes available. 

> : 

^ \ There has been considerable recent interest in the status of the first row CKM unitarity 

CN ' relation, \V U( i\ 2 + \ V US \ 2 + \ V u b\ 2 = 1. Recent determinations of \V US \ from (i) r^ 2 /r^ 2 , 

analyzed using lattice results for f K / f n as input [1,2] and (ii) data, using updated 
results for the kaon lifetimes, partial widths and form factor slope parameters, 
combined with lattice and ChPT work on the deviation of ff n (0) from 1 [3], yield results 
compatible, within errors, with 3-family unitarity. A major component of the error in 
both approaches is theoretical, and associated with uncertainties in the extrapolation of 
lattice data to physical values of the light quark masses. Here we discuss an alternate 
approach to extracting |V« S |, based on analyses of flavor-breaking sum rules involving 
hadronic X decay data [4, 5]. An important feature of this method is that the associated 
theoretical errors are completely independent of the chiral extrapolation uncertainties 
^ ■ which dominate the and Kya vs. llya. approaches. 
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The T-based method works as follows. With IlyL.^. the spin J parts of the flavor 

J V/A-ij ' 



ij = ud, us vector/axial vector correlators, p^,... the corresponding spectral functions, 



and R v / A;i j = T[x — > v T hadrons v / A;i j (y)]/T[r — > v r v e e (7)], x decay kinematics yield 
R v/A;ij = \2n 2 \V l} \ 2 S EW J\y t [w T (yr)p {0+l) (s) + w L (y T )p^(s)\ , (1) 

where y t = s/m 2 , Sew is a short-distance electroweak correction, wj{y) = (1 — y) 2 (l + 
2y), wl(v) = — 2v(l — y) 2 , and the superscript (0+1) denotes the sum of J = and 7=1 
contributions. As written, the RHS of Eq. (1) can be recast using the general FESR 
relation, / ?0 dsw(s)p(s) = ^ §^ =S(j dsw(s)Tl(s), valid for any fl without kinematic 
singularities and any analytic w(s). Rescaling the experimental decay distribution in 
Eq. (1) by (1 — Vt)V™ (k,m > 0) before integration leads to the generalized spectral 

integrals, Ry'^p and corresponding "(k,m) spectral weight sum rules". Analogous 
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spectral integrals and FESR's can also be written down for so < m\, for the separate 
correlator combinations Uy^j(s), sUyj^^s), and for general w(s). The corresponding 
spectral integrals are denoted generically Rfj(so) below. At sufficiently large sq the OPE 
representation can be employed for Tl(s) on the RHS of the FESR relation. 

With this notation, choosing a particular so, w(s) and correlator combination, the 
spectral integral difference 8R w (s ) = [R w ud {so) /\V ud \ 2 ] - [RZ(s )/\V us \ 2 ] vanishes in 
the SU(3) flavor limit. Its OPE representation, 8Rq PE (sq), thus begins at dimension 
D = 2. Using this representation, and solving for \V US \, one has [4] 

\Vus\ = y / ^SoV [(R: d (so)/\V ud \ 2 ) - 8R w OPE (s )] . (2) 

An important practical observation is that, for a range of weights studied in the literature, 
8Rq PE (sq) is at the few to several % level of R™ d {so) /\V u d\ 2 ■ High-precision determina- 
tions of \V US \ are thus possible from modest-precision determinations of 8Rq PE (so), pro- 
vided errors on the ud and us spectral integrals are kept under control [4]. The following 
technical points are relevant to constructing reliable versions of this analysis: 

• Because (i) the integrated D = 2 OPE for the purely 7 = contributions displays 
extremely poor convergence for all sq accessible in hadronic x decay [6], and (ii) all 
truncation schemes considered in the literature produce strong violations of spectral 
positivity [7], it is necessary to work with FESR's based on the difference, All, of 
the ud and us J = + 1 correlators. Fortunately, for a combination of chiral and 
kinematic reasons, the 7 = contributions to the experimental spectral distribution 
are strongly dominated by the well-known K and % pole terms, and hence can be 
straightforwardly subtracted. The desired 7 = 0+1 spectral difference is then easily 
determined. (Small residual corrections are constrained by analyses of the strange 
scalar and pseudoscalar sectors [8].) 

• The series in a s for [AIT] ^2 is slowly converging at the spacelike point on \s\ = sq, 
for all sq < m 2 [9, 10]. This necessitates using weights which emphasize other 
regions of the contour, where | a s (Q ) | is smaller and the convergence of the corre- 
lator series improved. The (k, 0) spectral weights are disfavored from this point of 
view since the modulus of (1 — y) k+1 peaks more and more strongly in the space- 
like direction (y = — 1) as k is increased. The convergence of the integrated D = 2 
OPE series is already very poor for the (0,0) case [9, 11]. Fortunately, alternate 
weights, w\o, vt>io and W20 [11], displaying significantly improved convergence be- 
havior [5, 11], exist in the literature. 

• Lack of empirical information on higher D VEV's necessitates neglect of D > 6 or 
8 OPE contributions. This neglect can be tested for self-consistency by employing 
weights w(y), with y = s/sq, and checking output \V US \ values for any unphysical 
so-dependence, since integrated D = 2k contributions then scale as 1/sq relative to 
the leading integrated D = 2 term. Such unphysical so -dependence can also result 
from truncation of a slowly converging integrated D = 2 series. Testing the output 
for so-stability is thus very important, especially when, as for the (0, 0) spectral 
weight, convergence is slow and conventional truncation error estimates, as a result, 
are potentially suspect. To make such SQ-stability tests possible, it is essential that 
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experimentalists provide full covariance matrices, and not just quote correlated 
errors for a limited number of weights and for sq = m 2 only. 

With currently available CLEO, ALEPH and OPAL data [12], ud and us spectral 
integral errors (~ 1/% and ~ 3 — 4%, respectively) lead to ~ 1/4% and ~ 1.5 — 2% 
errors on \V US \. We have employed the ALEPH data and covariance matrices, with the 
distribution for each exclusive mode rescaled by the ratio of the current world average to 
original ALEPH branching fraction for that mode [13]. The largest rescaling occurs for 
the K^K^ll^ mode, where the OPAL and CLEO results are in good agreement, but in 
disagreement with ALEPH. We quote results below for both the 3-fold ALEPH-CLEO- 
OPAL ('ACO') and 2-fold CLEO-OPAL ('CO') averages. The dominant us spectral 
integral errors will be dramatically reduced once analyses of the much larger data sets 
available from the B factory experiments become available. 

OPE errors are dominated by uncertainties in the leading integrated D = 2 series. The 
results reported below correspond to the PDG04 input, m s (2 GeV) = 105 ±25 MeV, 
whose square sets the overall scale for the D = 2 series. Other OPE input is detailed 
in Ref. [5]. The D = 2 error has two main sources: the truncation of the D = 2 series, 
and the overall [m s {2 GeV)] 2 scale uncertainty. The latter can be easily reduced through 
near-term improvements in our knowledge of m s , while the former are not subject to 
near-term improvement. A useful indicator for the quality of the convergence of the 
integrated D = 2 series is r k (so), the difference between the 0(af)-truncated correlator 
and Adler function versions of the D = 2 integral, scaled by the corresponding correlator 
integral. The two expressions must be identical to all orders but, at fixed order k, differ by 
terms of 0(a^ +l ). For a well-converged series, r^so) should be small, and decreasing in 
magnitude with increasing k. For wio, vt>io and w 2 o, we find ^^(m 2 ) = —0.06, —0.03, 
-0.01; r 2 3 A (m 2 t ) = -0.07, -0.05, -0.03; and r 2 3A{m\) = -0.08, -0.05, -0.03, re- 
spectively, confirming the (by design) improved convergence of the corresponding inte- 
grated D = 2 series. For the (0,0) spectral weight, in contrast, ^^(m 2 ) = 0.06, 0.20, 
0.67, confirming the poor convergence of the integrated D = 2 series. The latter result 
raises warning flags about the reliability of existing (0, 0) spectral weight truncation er- 
ror estimates (the difference between the correlator and Adler function evaluations alone 
is a factor of ~ 2.5 times larger than the original estimate of Refs. [4]). D = 2 errors turn 
out to be dominated by the (sizeable) truncation uncertainty (and hence not improvable) 
for the (0, 0) spectral weight, but by the overall scale uncertainty (and hence easily im- 
provable) for w\q, vt>io and w 2 q [5]. A more detailed breakdown and discussion of OPE 
errors can be found in Ref. [5]. 

The results of the so-stability tests for the various weights are displayed in Table 1 . 
One sees poor stability, and hence significant additional theoretical uncertainty (much 
greater than the total ±0.0009 theoretical uncertaintiy estimated in Refs. [4]), for the 
(0, 0) spectral weight, but good stability for w\q, m>io and w 2 q. 

Our final results are based on the w\q and w 2 q weights, and sq = m 2 , for reasons 
discussed in Ref. [5]. We find, for the ACO (CO) treatments of the us data, 

|K, I = 0.2223 (0.2246) ±0.0032 exp ±0.0038^ , (3) 

both results in agreement, within errors, with those of other methods. We stress that 
(i) the errors here are entirely independent of those of the other methods; 
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TABLE 1. \V US \ as a function of so for various FESR weight choices and the ACO and 
CO treatments of the us data, sq is given in units of GeV 2 



so 


w ACO 


<° 


<° 


<° 


w co 


W 10 


W 20 




2.35 


0.2149 


0.2220 


0.2243 


0.2201 


0.2172 


0.2236 


0.2255 


0.2218 


2.55 


0.2167 


0.2218 


0.2235 


0.2203 


0.2192 


0.2236 


0.2250 


0.2223 


2.75 


0.2181 


0.2218 


0.2230 


0.2207 


0.2207 


0.2239 


0.2246 


0.2229 


2.95 


0.2193 


0.2220 


0.2227 


0.2211 


0.2219 


0.2243 


0.2245 


0.2235 


3.15 


0.2202 


0.2223 


0.2226 


0.2216 


0.2228 


0.2246 


0.2246 


0.2241 



(ii) the experimental errors will be subject to significant near-term improvement once 
data from the B factory experiments becomes available; 

(iii) the theoretical errors for the (0,0) spectral weight approach are dominated by 
D = 2 convergence problems, and not subject to future improvement, making this 
version of the analysis unfavorable; and 

(iv) for the alternate weights, the theoretical errors are dominated by the D = 2 
scale uncertainty and reducible to the ±0.0013 level (or below) once the uncertainty 
in m s (2 GeV) is reduced to ±5 MeV [5]. Sub-1% determinations of \V US \ from hadronic 
T decay analyses are thus a realistic goal, with a further reduction in errors possible 
through averaging with the output from other methods. 
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